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ABSTRACT 

Comparison of linewidths of spectral line profiles of ions and neutral molecules have been recently 
used to estimate the strength of the magnetic field in turbulent star-forming regions. However, the ion 
(HCO + ) and neutral (HCN) species used in such studies may not be necessarily co-evolving at every 
scale and density and may thus not trace the same regions. Here, we use coupled chemical/dynamical 
models of evolving prestellar molecular cloud cores including non-equilibrium chemistry, with and 
without magnetic fields, to study the spatial distribution of HCO + and HCN, which have been used 
in observations of spectral linewidth differences to date. In addition, we seek new ion- neutral pairs that 
are good candidates for such observations because they have similar evolution and are approximately 
co-spatial in our models. We identify three such good candidate pairs: HCO + /NO, HCO + /CO, and 
NO+/NO. 

Subject headings: ISM: molecules - ISM: clouds - ISM: magnetic fields - magnetohydrodynamics 
(MHD) - stars: formation - ISM: abundances 



1. INTRODUCTION 

Observations of large velocity dispersions in the line 
profiles of molecules (line widths wider than expected 
from thermal motions) in interstellar molecular clouds 
is a signature of localized supersonic tu rbulence, while 
obser vations of the Zeeman effect (e.g., ICrutcher et al.l 
19991) and the polarizat ion of dust thermal emission (e.g. 



Hildebrand et al.l I2009D indicate that these clouds are 



permeated by magnetic fields. Thus interstellar molecu- 
lar clouds need to be studied in the framework of mag- 
netohydrodynamics. 

For typical molecular cloud conditions (e.g., B ~ 
10/iG, density ~ 10 3 cm~ 3 , and ionized fraction < 10~ 6 ) 
magnetic forces, which are experienced directly by ion- 
ized particles, are transmitted to the bulk of the neu- 
tral gas molecules via frequent collisions. However the 
coupling of the magnetic field with the gas is not per- 
fect. Motions of the bulk neutral gas with frequen- 
cies higher than the typical ion-neutral collision fre- 
quency can not be transfered to the ions. This rel- 
ative motion between cha rged and neutral part icles, 
known as ambipolar drift ()Mestel fc Spitzerl I1956D . re- 
sults in the dissipation of hydromagnetic waves in 
small scales flB raginsk iH 119651 : iKulsrud fc Pearcel 119691 : 
iZweibel fc Josafatssonl Il983[ ) and can cause differences 
between the velocity dispersion s pectra of the ions a nd 
the neutrals, as was recog nized bvlHoude et al.l (|2000al lbT) 
and explored further by iLi fe Houdd (|2008lV Indeed, 
iMouschovias et ail (|2011l) . in a comprehensive study of 
magnctohydrodynamic waves in weakly-ionized gases, 
showed that most hydromagnetic waves with wavelengths 
below a critical value cannot propagate in the neutrals 
because they are damped rapidly by ambipolar diffusion, 
but certain long-lived modes in the neutrals exist at all 
wavelengths; in contrast, waves in the ions damp very 
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rapidly at short wavelengths. 

Observationally, the narrower line profiles of ions com- 
pared to that of co-exist ing neutral species w as first no- 
ticed and investigated bv lHoude et al.l (|2000alfbl l . who ar- 
gued that the reason of this ion line narrowing effect was 
due to the presence of magnetic fields in turbulent molec- 
ular clouds that are on averag e not aligned wit h the local 
flows. To test this analysis, lLai et ail (|2003f) observed 
DR21(OH) at a high spatial resolution in H 13 CO+ and 
H 13 CN and confirmed that the ion-neutral pair is co- 
existent within their observing scale (6") and that the 
ion line widths are indeed n arrower than the neutral line 
widths. ILi fc Houdd ()2008l ) re-visited this magnetically 
induced ion line narrowing effect. They mapped the Ml 7 
star forming region in HCO+ (J = 4 -> 3) and HCN 
(J = 4 — > 3) and fitted the lower envelopes of the ve- 
locity dispersion spectra of the observed line profiles of 
the ion and neutral pair to a Kolmogorov-type power- 
law. They determined the ambipolar diffusion length 
scale using the fitting parameters, and used it to calcu- 
late the strengt h of the plane - of-sky component of the 
magnetic fie l d. IHezareh et al.l (|2010f) further tested the 
ILi fc Houdd ((20080 technique by mapping DR21(OH) in 
the optically thin H 13 CO+ and H 13 CN ( J = 4 -> 3) lines 
and obtained the turbulent ambipolar d iffusion length 
and m agnetic field strength in that source. IHezareh et al.l 
(2012) used interferometry maps of the ground transition 
of the latter pair of molecules in a few massive dense cores 
in the Cygnus X region to verify that the method is re- 
producible for both single-dish and interferometry data. 

One concern in interpreting these observations is that 
chemical differentiation might also be responsible for dif- 
ferences in the spectral appearance of different molec- 
ular species: species with different spatial distribu- 
tions would also exhibit different velocity profiles. Al- 
though HCN and HCO + have been shown to be approx- 
imately co-existen t in the observed m ol ecular clouds on 
the scales probed (|Houde et al.ll20d0aH b1: lLai et all 120031 : 
IHoude et al.l 120041 ) . their chemical co-evolution has yet 
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Fig. 1. — Line types and colors used to denote each of the mod- 
els studied. Solid normal-thickness red line: "reference" magnetic 
model; solid normal-thickness blue line: "reference" non-magnetic 
model. Dotted lines: "fast" models; dashed lines: "slow" models. 
Brown/purple lines: magnetic/nonmagnetic models with temper- 
atures differing from the "reference" models. Orange/cyan shaded 
areas: variation in C/O ratio. Thin/thick solid red/blue lines: 
lower/higher cosmic-ray ionization rate magnetic/non-magnetic 
models (see text for details). 

to be examined from a theoretical perspective. Here, we 
use both magnetic and non-magnetic dynamical models 
of core collapse coupled with, and self-consistently follow- 
ing, non-e quilibrium chemist ry both in the gas phase and 
on grains ([Tassis et al.H2012l hereafter paper I) to exam- 
ine whether the typically observed neutrals and ions are 
indeed theoretically expected to be similarly distributed 
in a star-forming core. We also repeat this exercise for 
other frequently observed neutral and ionized molecules, 
and we propose neutral-ion pairs that, according to our 
models, best maintain a co-spatial distribution. 

2. MODELS 

We are interested in the co-spatiality of various observ- 
able neutral-ion pairs as a function of density in a star- 
forming region. Since the behavior of the abundance of 
many molecules with density is sensitive to various model 
parameters, we use all Paper I models to assess the be- 
havior of the abundance ratios of neutral-ion pairs with 
density. The best candidate molecule pairs for studies of 
systematic differences of the linewidths between neutrals 
and ions will be those with abundance ratios that vary 
weakly with density for any set of model parameters. 

The physics of the chemodynamical models we use and 
its effect on molecular abundances are exhaustively dis- 
cussed in Paper I. Here, we briefly review the parameters 
we have varied for each class of models: temperature, 
C/O ratio, cosmic ray ionization rate, and a parame- 
ter controlling the time available for chemical evolution. 
The latter is the mass-to-magnetic-flux ratio for magnetic 
models, and the collapse delay time for non-magnetic 
models (an initial time period during which chemistry 
evolves but the core does not evolve dynamically, repre- 
senting an early stage of support due to turbulence which 
later decays). 

The initial parameter values for our "reference" mag- 
netic model are a mass-to-mag netic-flux ratio equal to 
the c ritical value for collapse (M ouschovias fc Spitzerl 
1976), a temperature of 10 K, a C/O ratio of 0.4, and 
the canonical value for the cosmic ray ionization rate 
( = 1.3 x 10 — s — 1 . The relevant parameter values of 
our "reference" non-magnetic model are a collapse delay 
time (hereafter "delay" ) of 1 Myr, with the other param- 
eters being identical to the "reference" magnetic model. 



For magnetic models we vary the initial mass-to- 
magnetic-flux ratio, examining two additional values: 
1.3 times the critical value (a faster-evolving, magneti- 
cally supercritical model), and 0.7 of the critical value 
(a slower, magnetically subcritical model). For non- 
magnetic models we examine two additional values of 
delay: zero, and 10 Myr. 

For each of these six dynamical models, the carbon- 
to-oxygen ratio is varied from its "reference" value by 
keeping the abundance of C constant and changing that 
of O. Additional values of C/O ratio we examined are 1 
and 1.2. This way, we study a "matrix" of 18 combina- 
tions (9 magnetic and 9 non-magnetic) of available times 
for chemical evolution and C/O ratios. 

To test the effect of temperature (T), we have varied 
each of the six basic dynamical models by changing T 
by a factor of ~ 1.5 from its reference value of 10 K and 
examined models with T = 7 K and T = 15 K. This 
results in 12 additional models (6 magnetic and 6 non- 
magnetic) . 

To test the cosmic ray ionization variation, we stud- 
ied four additional models (two magnetic and two non- 
magnetic) , which have the "reference" value for the tem- 
perature, C/O ratio, and mass-to-flux ratio or delay (for 
magnetic and non-magnetic models respectively), but for 
which £ is varied by a factor of four above (C = 5.2 x 10~ 17 
s- 1 ) and below (C = 3.3 x 10~ 18 s" 1 ) its "reference" 
value, covering the r a nge of observat i onal e stimates (e.g., 
IMcCall et all pffll : iHezareh et all (pOO^ L 

Figure Q] shows a visual representation of these suite of 
models and is a quick reference guide for the line colors 
and types we have used to depict each model. 

3. RESULTS 
3.1. HCO+/HCN 

Figure [2] shows the results of our aforementioned mod- 
els for the evolution of the commonly observed ion- 
neutral pair HCO+/HCN. We show the evolution of the 
central abundance ratio with central density (left col- 
umn) , and radial profiles of the abundance ratio at two 
different evolutionary stages (of central density 10 4 and 
10 6 cm -3 , central and right columns, respectively). The 
top row shows the effect of dynamics on the abundance 
ratio, while the second, third, and bottom rows demon- 
strate the effect of the CR ionization rate, C/O ratio, 
and temperature, respectively. 

As shown in the left column, the central abundance ra- 
tio, especially in the magnetic models, is not particularly 
stable with central density. The reason, as can be seen 
in Figure 5 of Paper I, is that although both HCO + and 
HCN deplete with increasing central density, in our mag- 
netic models HCO + depletes initially faster than HCN 
and therefore their ratio declines by almost three orders 
of magnitude by a central density of 10 4 cm -3 . However, 
at higher central densities, the depletion rate of HCO + 
slows down while HCN keeps depleting at its initial rate. 
This effect is responsible for the increase of the abun- 
dance ratio at higher densities. In non-magnetic models, 
the dependence of the abundance ratio on central den- 
sity is much less pronounced, as the behavior of each 
of HCO + and HCN with central density is qualitatively 
similar (Figure 5 of Paper I) , with only a small rise in the 
abundance ratio above 10 6 cm -3 due to some flattening 
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Fig. 2. — Left column: Evolution of central abundance ratio of HCO+ and HCN as a function of central number density. Middle and right 
columns: radial profiles of the HCO" 1 " / HCN abundance ratio, for two different evolutionary stages as quantified by the central number 
density. Line types and colors are as described in Fig. [T] Top row: Different dynamical models for fiducial values of the C/O ratio, CR 
ionization rate, and temperature. Seond row: effect of the CR ionization rate for the fiducial magnetic and non-magnetic models. Third 
row: effect of the C/O ratio for each dynamical model, for fiducial values of the CR ionization rate and temperature. Bottom row: effect 
of temperature for each dynamical model for fiducial values of the CR ionization rate and C/O ratio. 
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of the depletion rate of HCO + at the highest densities 
examined. 

This behavior is also reflected in the radial profiles of 
this abundance ratio, as seen in the middle and right 
columns of Figure [2] The magnetic models and the slow 
non-magnetic models (blue, cyan, and purple dashed 
curves) exhibit a remarkable variation between the outer 
layers of the star-forming core where significant amounts 
of HCO + still exist, and the innermost and densest parts. 
The fast non-magnetic models, however, have a much 
milder dependence on radius. At a more advanced evo- 
lutionary stage (n c ~ 10 6 cm -3 ), the minimum seen in 
the evolution of the central abundance ratio with den- 
sity in magnetic models also appears in the radial profile 
at intermediate scales. At scales comparable to the size 
of the core, the radial profiles of the HCO+/HCN abun- 
dance ratio flatten out for all models. 

The total variation of the abundance ratio across mod- 
els, parameters, and densities is large (several orders of 
magnitude). As seen in the top row of Fig. [2l most of 
the total variation has its origin in the dynamics, with 
the CR ionization rate having the second most impor- 
tant contribution (second row). However, different mod- 
els represent different possible clouds in nature, and for 
this reason the variation of the abundance ratio across 
models is not the relevant concern for the Houde et al. 
type of analysis; instead, the quantity of interest is the 
total variation of the ratio within a single model, across 
space and time. Still, as seen in for example in the left 
column of Fig. [21 the maximum possible variation within 
a single model can be as large as several orders of mag- 
nitude. 

Observationally, the intensity maps of the 
HCO + /HCN pair show highest degree of correla- 
tion at large length scales, about 0.5 — 1.0 pc (|Lo et al.l 
2009), which correspond to scales comparable with or 
larger than the entire core in our models (r/Ro > 1 in 
Figured]). However, this correlation can be reduced at 
smaller scales of about 0.1 pc, i .e., resolutions achie ved 
by interferometric observations (jCsengeri etaDHoiil). 

HCO + and HCN follow independent chemical networks 
and thus have no reason to be co-evolving. The two 
molecules evolve differently and cannot be considered co- 
spatial in most of our models. HCO + is more abundant 
at lower densities and larger spatial scales. Furthermore, 
the critical densities of the commonly observed rotational 
transitions of HCO + are lower than these of HCN. Con- 
sequently, HCO + would in principle be expected to trace 
larger spatial scales, regions with higher levels of turbu- 
lence, which would lead to wider linewidths for HCO + 
rather than HCN. However, thi s is contradictory to ob- 
servations. iHoude et al.l (|2000af ) observed the J = 3 — » 2 
and J = 4 — » 3 transitions of HCO + and HCN in sev- 
eral molecular clouds where they found these species co- 
existent, with the ion lines indeed more opaque but nar- 
rower than the neutral lines. This can be explained by 
the effect of the magnetic field on the dynamics of molec- 
ular ions through the cyclotron interaction, which can 
in turn be interpreted as a signature of ambipolar dif- 
fusion (jHoude et al.l 12004] ) and be used to estimate the 
strength of the magnetic field. Indeed, it was later shown 
that clouds with stronger magnetic fields (projected in 
the plane of the sky) show larger differ ences between 
the ion and neutral spectral linewidths f)Li et al.l 120101 : 



iHezareh et al.ll2012T) . In contrast, chemical and radiative 
transfer effects act in the opposite direction causing the 
widths of neutral spectral lines to become narrower than 
those of ions, thus underestimating the strength of the 
magnetic field. 

To make a more robust comparison of our models with 
existing observations we must perform radiative transfer 
calculations for the lines observed. We will return to this 
problem in the future. An observational correlation sim- 
ilar to HCO+ and HCN is observed for the H 13 CO+ and 
H 13 CN pair but since our chemical models do not incor- 
porate the chemistry of the isotopes of carbon and other 
elements, we do not include these isotopologue pairs in 
our discussion. For the remainder of this paper, we iden- 
tify additional ion and neutral species that are good can- 
didates for being considered as co-existent pairs. 

3.2. HCO+/NO 

Figure |3] shows the same pattern as in Figure^ i.e., the 
evolution of central abundance ratio with central density 
in the left column, and radial profiles of the abundance 
ratio for central densities of 10 4 and 10 6 cm -3 in the mid- 
dle and right columns, respectively, for the HCO + / NO 
pair. The central abundance ratio of this pair is signifi- 
cantly more stable with central density within any single 
model than the HCO+/HCN ratio, as HCO+ and NO 
evolve in a similar manner with central density, with the 
exception of the high C/O ratio models (third row: yel- 
low and cyan lines for the magnetic and non-magnetic 
models, respectively). In all other cases, the variation 
of the abundance ratio with central density within a sin- 
gle model is < an order of magnitude over the central 
density range we examine. 

A similar behavior is reflected in the radial profiles of 
the HCO + / NO abundance ratio, with most models ex- 
hibiting little variation. The large variation in the HCO + 
/ NO pair for the magnetic high C/O-ratio model (yel- 
low curve) only affects very small scales, which are also 
associated with high densities. At these densities, how- 
ever, the depletion of these molecules is significant (Fig. 
5 in Paper I), so these regions are not expected to have 
a considerable contribution to the observed lines. 

3.3. HCO+/CO 

The plots for the HCO + / CO pair is shown in Fig- 
ure 131 This pair exhibits very small variations of the 
central abundance ratio with central density for all but 
the fastest non-magnetic models. However, even in the 
case of these fast models (dotted blue, cyan, and purple 
curves), the variation is less than an order of magnitude 
for central densities above 10 4 cm" 3 . 

In the radial profiles, some variation in the abundance 
ratio occurs at large radii, but it is less than about an 
order of magnitude, at least in early times, and not much 
more than an order of magnitude in later times, even in 
the most-varying, fast non-magnetic models. 

It is not surprising that the abundances of the two 
molecules are co-evolving, since the dominant reaction 
for the production of HCO + in our models involves CO: 

CO + H+ ^HCO++H 2 , (1) 
while the destruction of HCO + proceeds primarily 



through dissociative recombination with an electron: 

HCO + + e- -» CO + H, (2) 

which leads back to creating CO, forming a produc- 
tion/destruction loop for the two molecules. 

3.4. NO+/NO 

Figure [S] shows the same plots as before for the NO + 
/ NO pair. This is the best example of co-evolution 
and chemical co-spatiality of an ion-neutral pair that we 
have identified in our coupled chemical/dynamical mod- 
els. The variation of the central abundance ratio with 
central density in any single model is at most a factor 
of a few. These findings are also reflected in the radial 
profiles as well. 

The reason for the co-evolution of the two molecules is 
that NO + is mainly produced through ionization of NO 
or charge- exchange reactions, also involving NO: 

NO + H+ ->• NO+ + H (3) 

and 

NO + C+ -> NO+ + C . (4) 

The ion is destroyed through dissociative recombina- 
tion with an electron: 

NO + +c~^0 + N, (5) 

so the abundance of NO + is coupled to that of NO. 

4. DISCUSSION 

In addition to the co-spatiality of the three new ion- 
neutral pairs that we have identified, their observabil- 
ity and critical densities for excitation of their rotational 
transitions must also be considered to decide whether 
they are viable candidates for studies of the ambipolar 
ion-neutral drift in star-forming regions. 

The most frequently observed transitions of HCO + 
have critical densities^ of the order of 10 5-6 cm~ 3 . NO 
generally has much lower critical densities, so overall 
unless higher critical density transitions of NO can be 
observed, this would pose a problem in the use of the 
HCO+/NO pair. The reason is that NO would trace 
lower-density, higher-turbulence medium, which would 
lead to wider linewidths for NO, a trend that could mimic 
the effect observed in the HCO + /HCN pair. Possible 
higher critical-density transitions would, for example, be 
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the J = 7/2 — > 5/2 transition of NO at 350.6 GHz, with 
a critical density of ~ 10 5 cm~ 3 . The multiplets of the 
first two rotational transitions of NO, i.e., J = 3/2 — > 1/2 
and J = 5/2 — > 3/2 have been obse rved in both dark and 
massive clouds (|Gerin eFaT][T99l at 150 GHz and 250 
GHz with relative abundances calculated to vary from 

io- 7 -io~ 8 . 

A similar problem arises with the HCO+/CO pair, as 
CO, with generally low critical densities in usually ob- 
served transitions, traces more extended regions. Be- 
cause the abundance of CO is high, such transitions 
would result in saturated and optically thick lines. To 
avoid this effect, one needs to use higher CO transitions. 
A possibility would be the J = 6 — >• 5 transition at 691 
GHz, with critical density around 10 5 cm' 3 . 

The NO + /NO pair is in principle the best choice 
in terms of co-spatiality of the neutral and the ion, 
but although NO + has sub-milli meter and radio trans i- 
tions identified in the laboratory (jBowman et~aII[T98^ . 
the ion has not yet been observed in molecular clouds. 
This may be due to a general trend in which nitrogen- 
bearing molecules, with the exception of NH3 and N2H + , 
have not been targeted for observations as much as the 
carbon-bearing species, in part because nitrogen chem- 
istry in molecular clouds is even less understood (espe- 
cially at low tempera tures) than that of carbon (e.g., 
iHilv-Blant eTaI|[20Toh . 

We note that our models do not include the effect of 
turbulence; as a result, the relative abundance variation 
that we find is an upper li mit, since turbu lent mixing 
would lessen this variation ()Xie et al.lll995h . Also, our 
models represent prestellar cores, and do not include the 
effect of feedback sources, because our parameter studies 
do not feature very large variations of the temperature, 
and UV ionization has not been accounted for. However, 
regions with feedback sources are avoided when applying 
the Houde et al. method. 
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